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Abstract 

We study the rare kaon decays in the framework of general SUSY models. Unlike the results 
in the literature, we find the contributions from the gluino exchange to the branching ratio 
of K + — > 7i + uu can reach the central value (~ 1.5 x 1CU 10 ) of the new E787 data while the 
predicted value of standard model is less than 1CU 10 . We also find that the same effects also 
enhance the decays of K L — > 7r°z/z/, K l — > ir°e + e~~ and K L — > 



One of the essential reasons for the success of the standard model (SM) is that it naturally 
satisfies all measured phenomena of flavor changing neutral current (FCNC) in K and B 
meson systems. The FCNC processes are forbidden at the tree level and suppressed at the 
loop level by the Glashow-Iliopoulous-Maiani (GIM) mechanism and by the small quark 
mixing matrix elements which involve the transitions between the third and the first two 
generations in the Cabibbo-Kobayashi-Maskawa (CKM ) matrix [|IJ. Therefore, those models 
which have no GIM mechanism and flavor mixing suppression could largely enhance the 
FCNC decays. Through these decays, one can search for the existence of physics beyond the 
SM. 

Recently, the processes associated with \AB\ = 1 and \AK\ = 1 have further progresses in 
experiments. The representative of the former is the decay modes of B — > K£ + £~ {£, = e, /i), 
observed at the Belle detector in the KEKB e + e~ storage ring with the branching ratio 
(BR) of Br(B K£+£~) = (0.75igjf ± 0.09) x lO" 6 , while the SM expectation is around 
0.5 x 10~ 6 ||. The latter is K + n + uu. The new data on BR from E787 f§ at the 
Alternating Gradient Synchrotron (AGS) of Brookhaven National Laboratory is given by 



Br{K + -> tt+z/z/) = (1.57±J;g) x lO" 10 , 



but the predicted BR in the SM is (0.72 ± 0.21) x lO" 10 [§, §. 

By comparing the theoretical uncertainties in both processes, one can find that the BRs 
of B — > K£ + £~ involve three independent form factors and their errors can not be reduced 
by the exiting measurements. Therefore, before studying the new physics effects on such 
processes, it is necessary to know how large the true theoretical uncertainties are. On the 
contrary, K + — > 7i + uu has less errors from the hadronic matrix element and is free of the 
long-distance uncertainty 0. For eliminating the effect of K + — > tt + matrix element, the BR 
of K + — > 7i + uu can be related to that of K + — > 7r°e + z/ with measured BR of 0.0482 by using 
the isospin symmetry. Even the corrections to the isospin limit have also been included by 
Ref. ||. Besides K + — > ix+vv, the relevant decays, such as K L — > fi + fi~, K L — > 7r°i/z/, and 
Kl — * 7r°e + e~, have the similar characters; especially the BRs of the last two are related 
to CP violation (CPV). By means of loop induced effects and less hadronic uncertainties, 
such kind of rare kaon decays provides good candidates to test the SM ||. On the other 
hand, it was pointed out by Ref. [|Kj that (3 or 1; one of the three angles in the CKM 
matrix, can be described by combining the BRs of K + — > n + uu and —>■ 7r°z/z/, denoted by 
[sin 20i] 7W p. As well known, this angle is also related to the time-dependent CP asymmetry of 
the B — > J/^/K s decay, expressed by [sin 20i] j/^k s - It is obvious [sin 2</>i] 7n , i j = [sin 20i] j/^k s 
in the SM. However, once introducing new CP phases, the identity will be changed so that 
[sin 2(0! + 6k)]ttvd 7^ [sin 2(0! + 0b)]j/vk 3 where 0k and 0b are the effects of new physics on 
K and B decays, respectively. Hence, by comparing sin 20! measured from rare kaon decays 
and the asymmetry in the B meson system, we can also tell whether there exist new physics. 

Super symmetric (SUSY) theory not only supplies an elegant mechanism for the breaking 
of the electroweak symmetry and a solution to the hierarchy problem, but possesses abundant 
flavor and CP structure. Besides the original CKM matrix, the SUSY models introduce new 
flavor mixing effects, such as the upper and down type squark mixing matrices. The new 
CP violating phases in SUSY models can arise from the trilinear and bilinear SUSY soft 
breaking A and B terms, the [i parameter for the scalar mixing as well as gaugino masses. 
Unfortunately, it has been shown that with the universal assumption on the soft breaking 
parameters, these phases are severely bounded by electric dipole moments (EDMs) |lTJ so 
that the contributions to e and e' are far below the experimental values. In the literature, 
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some strategies to escape the constraints of EDMs have been suggested. They are mainly 
(a) by setting the squark masses of the first two generations to be as heavy as few TeV 
fl2|l but allowing the third one to be light; (b) by including all possible contributions to 
EDMs such that somewhat cancellations occur in some allowed parameter space [T3J ]; 
and (c) with the non-universal soft A terms instead of universal ones. In particular, those 
models with non-universal parameters have been demonstrated that they can be realized in 
some string-inspired models J 
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In this paper, we 
> 7T + vv and other 



Moreover, without the universal assumption, 
the corresponding off-diagonal terms for the left-right mixing of the squark mass matrix 
are unnecessary to be proportional to the light quark mass directly [IE| 
will show the implication of the generalized A terms on the BRs of K + 
relevant rare kaon decays. 

For simplicity, we adopt the mass-insertion approximation in which the masses of squark 
are taken as degenerate approximately We will only concentrate on the Z^d^^P^^s 
effective interactions introduced by Z-penguin diagrams with Pl{r) — (1 T 7s)/2- Although 
box diagrams also contribute to s — > dvi>, one can easily check that the effects compared to 
those of the SM are suppressed by M w jM\ with M x being the typical SUSY mass scale. 
Here, we will not discuss them. As to the dipole vertices such as Z fJ, da^q u s, because they 
are suppressed by q/Mz and by the light quark masses m s (d), their contributions are also 
negligible. Hence, there are two mechanisms to generate effective couplings for Z^dj^Pi^s. 
One is from the penguin of chargino exchange [^U], E_T|, [22|, and another one is from the 
same diagram but with gluino exchange [EDI. Comparing to the contributions of chargino 



and gluino, due to smaller couplings, the effects of neutralino are always negligible. 

In the literature it is claimed that due to the magnitude of the relevant mass-insertion 



parameters being proportional to the light quark masses [p_l|, 22, 23 1 or proportional to the 
suppressed quark mixing matrix elements J2(J, the effects from gluino are much smaller 
than those of the SM. As a consequence, the dominant one comes from the chargino sector 
that the mass-insertion terms are associated with stop-quark effects and the corresponding 
constraints on the squark mixing matrix elements are looser |E2| . However, it is shown by 



the new data of E787 that the central value has been close to the predicted value of the 
SM. In spite of the still large errors, if the central value hints the existence of new effects, 
actually it implies that the models with very large contributions to K + — > 7i + uu should be 
further constrained. On the contrary, it increases the possibility for the new physics which 
is compatible with the SM. Due to the enormous progress in the SUSY flavor physics, what 
we want to emphasize is that the contributions from the Z-penguin with gluino exchange are 
not negligible and can have sizable effects if the conventional soft SUSY breaking A terms 
possess more general flavor structures or are non-universal (nondegenerate). To accomplish 
the purpose, in the following analysis, we consider the case that the contributions from 
chargino are small. 

We start by writing the effective interactions for Z coupling to quarks and squarks as 
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£ = -*—Y,[ftd»(T3 f -s 2 w Q f )f L 
cw - 



-r R d, (s 2 w Q f )f R ]z»-^- j:h 



c w 



(T 3/ - s 2 w Q f )P L - s 2 w Q f P R \fZ» 



where / denotes the fermion and /l(_r) corresponds to the superpatner of / with the chirality 
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L(R), s-wipw) = svn.6 w (cos9w) with 6 W being the Weinberg angle and Qf is the charge of 
/. On the other hand, the relevant squared down type squark mass matrix is described by 



M 



D 



[ m D>LR 



(3) 



= (M|) LL + m^-^(Ml + 2M^)i, 

2 \ f\/r2 



(m 6 ) LR = (M D ) LR -/x*tan/3m D , 
(m|) M = (M^ + m^-^^Mlsin 2 ^!, (4) 



where we have adopted the so-called super-CKM basis that the quarks have been the mass 
eigenstates so that is the diagonalized down quark mass matrix, 1 denotes the 3x3 unit 
matrix, the definition of the angle (3 is followed by tan/5 = v u /vd with v u and Vd being the 
vacuum expectation values (VEVs) of Higgs fields $ M and $ d responsible for the masses of 
upper and down type quarks, respectively, \i is the mixing effect of & u and & d , (M 2 ^) ll(rr) 
stand for the soft breaking masses for down type squarks and (M 2 ^)^ describe the trilinear 
soft breaking couplings and are written as 

(Ml) LR = l±V DL A*vi H , (5) 

where Vd l{r) transform the left(right)-handed quarks from weak eigenstates to mass eigen- 
states and Af- = Yf-Af- with Yf- and Af- being Yukawa and soft SUSY breaking matrix, 
respectively 

Compared to the My/ scale, due to the smallness of the involved momentum or momentum 
transfer and the masses of external legs, we drop them in our considerations so that one 
can easily get the conclusion by the similar situation to the gauge invariant requirement 
on 7 — s — d vertex that the one-loop contributions from the left-left (LL) mixing of the 
down squark mass matrix are vanished. Therefore, in order to get the effective interaction 
Z — s — d that the incoming and outgoing particles carry the same chirality, it needs a double 
mass-insertion in the squark propagator, i.e., the appearance of ABiM^^) B a with 



A(B) = L(R) or R(L) and j = 1, 2, 3 is necessary p0[ . In general, we also consider g — qi — QL 
and g — <1r — qR vertices simultaneously. The situation is different from the case of chargino 
exchange in which only the left-handed couplings give the main contributions. 

According to the interactions in Eqs. (§) and (||), by including the self-energy diagrams 
and all possible emissions of Z-boson from the propagators of internal squarks, the effective 
interactions for s — > dZ can be derived as 

C~ g = C z [Z L {x)d^P L s - Z R {x)d^P R s]Z» (6) 

with 

l^Z = — ^-rM 7 



V2tt 2 Z sin^ 
s sin 2 
I2a f 



Za(x) = as ^ 9w C F G{x) £ (4W^W (7) 



J=l,2,3 
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where Cp = 4/3, (6^)ab = (^I^Oab/w-, rriq is the average mass of squark in the super- 



CKM basis and 



G{x) 



2x 2 + 5x — 1 3x 2 \nx 



2(x - l) 3 



(x-1) 



with x = rri~/rri~ and rrig being the gluino mass. We note that according to Eq. (|6]) the 
effect of the different chirality is opposite in sign each other. As known, the associated 
matrix elements for relevant kaon decays is < Tr\d^f^s\K > so that the opposite sign actually 
reflects somewhat cancellation between different chiral couplings. We will see later that in 
some SUSY models, the cancellation is almost complete. 

Altogether, the effective interactions combined with those of the SM can be written as 



C 



(x S M(xt) + Z L (xfj d^^P L s 



■Z R (x) dj^Pns 



with 



X S M{x t ) = A C P + X t X (x t ). 



The explicit expressions of functions X (x t ) and Y (x t ) can be found in Ref. |23[. Pq is given 
in Ref. p|] . From the formulas in Ref. |23|] , the BR is expressed as 



Br{K J 



7T VlV X) 



{ImFf + {ReFf 



(9) 



with 



X~ g (x) 



jPo + ^Mxt) + y 5 X,(x) 
3a 2 Br(K+ -> 7r°e+z/) 8 

r K+ o 2 • 4/) ^ - 

2ir z sm Uw 
Z L {x) - Z R (x), 



(10) 



where x t = m 2 /m 2 v , Aj = V* s V ic i with A c being real to a very high accuracy, and r^M 
0.901 summarizes isospin breaking corrections in relating K + — > -k + vv to K + 
Here, the isospin breaking corrections from the quark mass effects and electroweak radiative 
corrections have been calculated in Ref. ||. By means of Eq. (|I~0D , we clearly see that if 
the SU(2)l breaking effects for left- and right-handed coupling are the same, the influence 
on Br(K + — > 7i + uu) also vanishes. With the same effective interactions, the contributions 
from the gluino exchange to BRs of Kl — > 7r°e + e~, Kl — > t^vv and Kl — > yU + /i~can also be 
described by 



Br(K L 



TT°e + e ) dir 



Br(K L -> 7r°m) 



Br{K L 



M A* )SD 



\Imy 



7A\ 



-K L {ImFf 



\Imy 



7V\ 



'ID 
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with 



and 



Via = ^Y (x t ) + ^Xg(x), 
Vtv = ^[l + (l-Asm 2 e w )C (x t )] 
+^(l-4sm 2 e w )Xg(x), 



= r £L 3 a ^r(^ 7 A^) A8 = Lg9 , 10 _ 10 



li 27r 2 sin 4 6^ 

where k = «l/6, = 1.68 x 10~ 9 , and C (x t ) and A c from the charmed loop can be found 
in Ref. @. 

The essential question is whether the involving gluino contributions can yield the BR of 
K + — > 7r + z/z/ as large as that given by E787. To explore the possibility, we have to analyze 
the constraints on the relevant parameters. For simplicity, we just consider the case for j=3 
in Eq. ( PH) so that the involving mass-insertion parameters are only (<5| 3 )ab and (S^ab 
although this assumption is unnecessary. To estimate the CP violating effects, we take that 
the real and imaginary parts of relevant parameters are approximately the same in order 
of magnitude, such as \Im(5fj)AB\ ~ |-Rc(^-)ab|/2. Under our assumption, it is known 
immediately that the bounds on the parameters are from BriB —>■ Xg'j) and the — B^ 
mixing. The former constrains (5$ 3 )ab while the latter is {S^ab- In addition, it is worth 
mentioning that charge and color breaking (CCB) minima and the potential unbounded from 
below (UFB) may also give strict bounds ||25|| . To relax the constraints from the vacuum 



instability, we adopt the following two strategies: (a) According to the result in Ref. [g5 



the model independent upper bounds on (S^)lr can be described by m&y2 + m~ /m~/mq. 
By taking the slepton mass to be few TeV, the constraints are compatible with those directly 
obtained by FCNC decays, (b) Our universe is resting a false vacuum. From the analysis 
in Ref. P^] , the constraints from the conditions of CCB minima and UFB can be relaxed if 
the life time of the metastable vacuum is as long as that of the present age of the universe. 



Although the requirements of CCB and UFB in [25] are necessary, after all, they are not 
sufficient. Hence, in our numerical calculations, we use the constraints gotten from FCNC 
processes. 



Since the relevant constraints have been considered in Ref. [27], we display the bounds in 
Table |Tj for fixing the specific chirality. The values in the entries of Table [I| are for rriq = 500 
GeV. For the different choices, the values for the second column need to be multiplied by 
mg/500 but it is (mg/500) 2 for the third column. The parameters that all R(L) are replaced 
by L(R) depend on the details of SUSY models. As illustrations, we show three possible 
situations in SUSY models: 

• Scenario I: Zl — Zr. This approximation is equivalent to (8fj)RL — (Sf^RL- It has been 
shown recently that if the Yukawa and soft-breaking A matrices are hermitian, the identity 
could be realized |28], 0. As a result of Eq. flTU|), the contributions to rare kaon decays 



all vanish. Nevertheless, this scenario implies that the hyperon CP asymmetry in SUSY is 
one order of magnitude larger than that of the SM ]3l] . 



5 



Table 1: The bounds on the relevant parameters [27 



x \Re(5^) RL \ x 10 2 


\MS^)lr\ x 10 2 


G(x) 


0.3 4.57 


1.16 


0.36 


1.0 3.81 


1.43 


1/4 


4.0 4.16 


2.68 


0.18 



• Scenario II: \Z R — Z L \ ~ \Z L ^\. Besides the negligible contributions to rare K decays 
in scenario I, another shortcoming is that the predicted e' is also far smaller than that of the 



experimental measurement p9j. In order to deal with the small e' problem and escape the 



constraint from the EDM, the asymmetric soft-breaking A matrix is proposed by Ref. |32 



and it is also found that such kind of the asymmetric property could be realized in some 
string-inspired supergravity models. Hence, the situation of scenario II can be reached if the 
SUSY soft-breaking A matrix is asymmetric. According to the results of Eqs. (§) and (|TT| ) 
and the constraints in Table [l], the BRs of relevant rare kaon decays are found as follows: 

Br(K + -> it + uu) = 1.55 x 10~ 10 , 

Br(K L -»■ tt°uu) = 2.50 x 10~ 10 , 

Br (K L -+ 7r°e+e-) dir = 3.47 x HT 11 , 

Br(K L -> (i + fi-)sD = 1-85 x KT 9 . (12) 

Because our purpose is to demonstrate that the gluino contributions from penguin diagrams 
to rare kaon decays could reach the current experimental ranges, the values in Eq. ([12]) are 
obtained by setting rriq = 640 GeV with x = 0.3 and the sign of Zum is the same as that 
of SM. And also we choose the value of X t such that the results of the SM for K + — > tt+vv, 
K L -> vr°z/z/, K L -> vr°e + e-, and K L -> are 0.58 x 10~ 10 , 0.46 x 10~ 10 , 0.7 x 10~ n and 

5.26 x 10~ 10 , respectively. 

• Scenario III: Z R <C Zl or vice versa. In this situation, as expected, the numerical 
predictions should be similar to the case of scenario II. 

Finally, we give a brief discussion on chargino contributions. According to the analysis 
of Ref. |[22|| , one can understand that unlike gluion case, the main effects of chargino on 
Zds effective interaction come from the wino component where only left-handed couplings 
are involved. That is, only (5%j) LR^ji) rl have the significant contributions. However, if 
the trilinear soft breaking Af, and Af- parameters arise from the same origin and have the 
same order of magnitude, due to the smaller weak couplings, the effects of chargino could 
be much smaller than those of gluino if both superparticles have the same masses. On the 
other hand, if we allow that A u is different from A d entirely, in order to guarantee that the 
effects of chargino are negligible, the mass of lightest chargino is taken as heavy as 0(TeV) 
so that although the constraints of relevant mass insertion parameters are not strict, all 
contributions from chargino will become insignificant. 

In summary, we have studied the rare kaon decays in the framework of general SUSY 
models. In terms of the rich flavor structure, we find that the effects of the gluino exchange 
can make Br(K + — > 7i + uu) up to the central value in the new BNL-E787 result. With the 
same effects, the BRs of the short distance contributions to Kl — > 7r°i/z/, Kl — > vr°e + e~ and 
K L — > fi + fi~ are larger than those in the SM. 
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